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We study the level statistics (second half moment 7o and rigidity A3) and the eigenfunctions of 
pseudointegrable systems with rough boundaries of different genus numbers g. We find that the levels 
form energy intervals with a characteristic behavior of the level statistics and the eigenfunctions in 
each interval. At low enough energies, the boundary roughness is not resolved and accordingly, the 
eigenfunctions are quite regular functions and the level statistics shows Poisson-like behavior. At 
higher energies, the level statistics of most systems moves from Poisson-like towards Wigner-like 
behavior with increasing g. Investigating the wavefunctions, we find many chaotic functions that 
can be described as a random superposition of regular wavefunctions. The amplitude distribution 
P{ip) of these chaotic functions was found to be Gaussian with the typical value of the localization 
volume Vioc ~ 0.33. For systems with periodic boundaries we find several additional energy regimes, 
where Io is relatively close to the Poisson-limit. In these regimes, the eigenfunctions are either 
regular or localized functions, where P(ip) is close to the distribution of a sine or cosine function in 
the first case and strongly peaked in the second case. Also an interesting intermediate case between 
chaotic and localized eigenfunctions appears. 
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I. INTRODUCTION 



^ I Quantum billiards are quite simple models for many practical applications in solid state and nuclear physics, as 

■ e.g. quantum dots, microdisk lasers and electron transport in microstructures. An important means for the study 
I of quantum billiards is the statistics of the quantum mechanical energy levels of a given system, where the distance 

■ distribution p(s) of the normalized distances s a = (E a+ i — E a )/(s) between two consecutive energy levels E a+ i 
I and E a with the mean distance (s) has the following two limiting cases, (i) The Poisson distribution, pp(s) = 

■ exp[— s], which is the distance distribution of uncorrelated numbers E a and (ii) the Wigner distribution pw( s ) — 
I 7rs/(2(s) 2 )exp[— 7rs 2 /(4(s) 2 )]. An example for the case (i) are the energy levels of a single symmetry group of an 

electron in a 2D potential well in the shape of a square, rectangle or circle (integrable systems). As a second example, 
localized states in disordered systems tend to the Poisson distribution with increasing system size, which is largely 
£h ■ used in solid state physics to distinguish between localized and extended states in disordered systems [1]. An example 
I for the case (ii) are the energy levels of chaotic billiards, as e.g. the stadium or the Sinai billiard. 

In close analogy is the behavior of a classical particle in a billiard, which undergoes elastic reflections at the walls. 
It also has two limiting cases, depending on the billiard geometry. If the billiard is chaotic, the motion of the particle 



jj^ ■ is ergodically extended over the whole energy surface in phase space. Two particles whose trajectories are very close 
at the beginning, diverge exponentially from each other. If the system is integrable on the other hand, the motion of 
the particle is restricted to a two-dimensional torus in phase space and neighboring trajectories diverge only linearly 
from each other. 

There are, however, several classes of intermediate systems between the two limiting cases, as e.g. polygonal 
pseudointegrable systems [2-4], on which we focuse in this paper, or systems with a pointlike scatterer ("Seba- 
Billiards") [5]. Like in integrable systems, the motion of a classical particle in a pseudointegrable system is restricted 
to a two-dimensional surface in phase space. However, these surfaces do not have the shapes of tori but are more 
complicated objects with more than one hole. They are called " multihandeled spheres" . Examples for pseudointegrable 
systems are polygons with only rational angles riiTr/mi, with n^m; £ N and at least one rii > 1. They are described 
by their genus number 

M 
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which is equal to the number of holes in the multihandeled sphere in phase space. Here, J is the number of angles 
and M is the least common multiple of the m*. The reason that those systems are not completely integrable is their 
property of beam splitting. At some points in their geometry, neighboring trajectories of particles can be split into 
two opposite directions (see Fig. 1). 

Several works [6-11] found the distance distribution of pseudointegrable systems intermediate between the Poisson- 
and the Wigner distribution. Ref. [6] presented a numerical investigation of pseudointegrable billiards of small genus 
numbers, where the boundary was approached by a step function, which was arranged along the curved shape of the 
(chaotic) Sinai billiard. With decreasing size and increasing number of the steps, the level statistics of this system 
approached the Wigner statistics. However, as it was pointed out in this work, the analysis was performed at very 
small energies, where the wavelengths were larger than the step sizes. So, it was assumed that the Wigner-like level 
statistics was not an inherent feature of the pseudointegrable shapes, but came from the underlying assymptotic shape 
of the Sinai billiard. Pseudointegrable systems of small genus numbers g — 2 and 3 were investigated numerically [7,8] 
and experimentally [8], showing an intermediate level statistics that clearly deviated from the Wigner distribution. 
In [9], numerical simulations on systems with increasing genus numbers up to g m 1000 indicated that for not too 
small energies, the distance distribution changes systematically from Poisson-like towards Wigner-like behavior with 
increasing g. Also for Seba billiards, it was shown by general arguments and numerical calculations that the level 
statistics is intermediate and close to Wigner-like behavior for small level distances [5,12-14]. The Wigner-like behavior 
is increased with the number and the coupling strength of the pointlike scatterers in the system [12,14]. 

The level statistics seems to be connected to the properties of the eigenfunctions. For example, in a billiard with 
g = 2 it was found that there are regular and irregular eigenfunctions coexisting [10]. In the Husimi representation, 
the eigenfunctions of systems with small g show signatures of pseudointegrability, whereas those of systems with large 
g tend to be irregular [11]. 

In this paper, we want to investigate the energy dependence of the level statistics of a special class of pseudointegrable 
systems with high genus numbers. At low energies, it is known that one can observe a level statistics that deviates 
from the high-energy limit, as shown in Ref. [6] for pseudointegrable and in Refs. [7,15-17] for chaotic billiards. Here 
we find that for our systems, there can be many energy windows, where the level statistics is comparatively close to 
Poisson statistics, and other energy intervals, where the behavior is close to Wigner statistics. We show that this 
behavior is correlated to the properties of the eigenfunctions and that there exist several characteristic types of such 
intervals, even at energy values, where the boundary roughness is resolved. 

The paper is organized as follows: In section II, we introduce the special pseudointegrable geometries that we 
consider. In section III, we show the results of the level statistics of the eigenvalues for several systems of different 
values of g. In section IV, we investigate the eigenfunctions in these energy intervals and introduce several 

quantitative measures, i.e., the localization volume the amplitude distribution P(ip^) and the behavior of the 

eigenfunctions in n, m-space. The functions ip and \& differ in their normalization, ip = ^/A^ , with the area A of the 
system. Finally, in section V, we apply those measures to many eigenfunctions in the different energy intervals. We 
find that in intervals where the energy-levels behave Poisson-like, the eigenfunctions are either localized or regular. 
Energy intervals with level statistics close to Wigner-like behavior, on the other hand, contain eigenfunctions that 
are random superpositions of plane waves. Also an interesting mixed case appears, where the superposition of plane 
waves leads to weak localization. 



II. SYSTEMS AND CALCULATIONS 



We consider a membrane of the sound velocity c that lies in the xy-pl&ne and vibrates in the z-direction. When the 
restoring forces are considered as scalar, the vibrations of this membrane are described by the Helmholtz equation 

A^ a \x,y) = -^ a \x,y), (2) 

with the a th eigenfunction ~$>( a \x, y), the corresponding eigenvalue us 2 a . The boundary conditions can be of Dirichlct 
or of Neumann type, referring to a membrane that is kept fixed at the boundary or that can vibrate freely, respectively. 
Equation (2) has the same form as the stationary Schrodinger equation with zero potential V = inside the system. 
Therefore, under Dirichlet boundary conditions, which refer to an infinite potential on the boundary, it also describes 
an electron of mass \i in an infinite potential well. In this case, one has to replace w^/c 2 by 2/j,E a /H 2 , with the energy 
eigenvalue E a . 
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For the numerical calculations, Eq. (2) is discretized on a square lattice, which reduces the problem to the diago- 
nalization of a symmetric matrix, which is carried out here by the Lanczos algorithm [18], a numerical procedure to 
compute eigenvalues und -vectors of sparse N x TV-matrices by reducing them iteratively to a tridiagonal form, for 
which effective algorithms exist. The eigenvalues E a are calculated numerically under Dirichlet and Neumann bound- 
ary conditions and their spectra are analyzed by the means of level statistics. The E a are dimensionless quantities, 
as we set h 2 / (2[i) = 1 and d = 1, where d is the lattice constant of the discrete lattice. 

As a model for our studies we chose rational billiards of the shape presented in Fig. 1 with different parameters, 
refering to different numbers, widths, distances and heights of the "teeth". The genus numbers of these geometries 
can be easily increased by increasing the number of teeth. Two different angles occur, ip — n/2 and ip = 3tt/2. 
Applying Eq. (1) we find that g = 1 + Gj, where Gi is the number of angles of values ip = 3ir/2. The systems have no 
symmetry axes and it is therefore not necessary to separate the calculated eigenvalues according to their symmetry 
groups (which would be technically difficult). The considered parameters range from systems with only few teeth and 
therefore very small genus numbers up to systems of g = 101. 



We now analyze the eigenvalues by the means of level statistics. The energy levels E a are normalized ("unfolded") 
to new values e a , such that their mean distance (s) is equal to 1. Then, we calculate the following two quantities: 
(i) From the nearest neighbor spacing distribution p(s), it has become common to calculate the second half moments 



which lie between the two limiting values /Signer ^ q 537 an j jPoisson _ ^ ^ gee j^ e f s [19^20]). This enables us to decide, 
if the statistics is closer to Wigner or closer to Poisson by comparing just one number I and is more comfortable in 
handling than p(s) itself. 

(ii) Another measure, which turns out to be even more sensitive to g is the spectral rigidity A 3 (L), where L is the 
length of the considered energy interval [21]. A 3 (L) starts from the integrated density of states N(e) = N(e a ) = 
^2n=i ®( £ a — £ n) of the unfolded energy levels e a , which is a staircase and can be approached by a straight line of 
slope one. A 3 (L) is defined as the least square deviation, 



where Min rii r 2 means that the parameters r\ and r2 are chosen such that the line n 4- T2£ is the best fit of N(e). 
For the calculation of A 3 (L) we use the technique derived in Ref. [22]. The limiting values are A 3 (L) = L/15 for 
intcgrable systems and A 3 (L) = ln(L)/7r 2 — 0.07/-7T 2 + 0(L _1 ) for the ensemble of Gaussian orthogonal matrices 
(GOE) [23], which serves as a general accepted good limit for chaotic systems. This means that in the first case, 
A 3 (L) increases linearly with L, and in the second case logarithmically. As discussed above, we expect intermediate 
behavior for pseudointegrable billiards. 

In this paper, we will basically use Iq as the easier of the two measures. However, we first want to compare the 
behavior of 1$ and A 3 for several systems, ranging from very small values of g — 3 until g ~ 49 at higher energies. 
The parameters of these systems are shown in Tab. I. All systems have roughly the same area and exactly the same 
widths and heights of the teeth, a = h = 8d, where d is the lattice constant. 

In Figs. 2(a) and (c) we plotted the values of Jo for Dirichlet and Neumann boundary conditions versus the energy. 
In Figs. 2(b) and (d) we show the A 3 (L)-data of the fixed energy intervals E{= uj 2 /c 2 ) G [1.5,2.0] versus L. The 
energy interval is chosen such that the wavelength, A = 2ttc/oj is smaller than the widths of the teeth and therefore 
the boundary roughness is resolved. The behavior of I shows many fluctuations but its average value decreases 
systematically with growing g towards the Wigner limit, which is shown as a solid line. Comparing Figs. 2(a) and 
(c), we also see that the behavior for both boundary conditions is quite similar. The only exception occurs at very 
small energies, where Iq lies considerably higher in the Dirichlet case. The reason for this is that the condition ip — 
exactly at the boundary prevents the long-wavelengths eigenfunctions to penetrate into the small boundary teeths. 
The same "screening" of the boundary roughness at small energies has already been observed for fractal drums [24]. 
For A 3 of pseudointegrable systems it was shown in [25] by a semiclassical periodic orbit theory that it depends on 
several details, e.g. on the energy interval and on the area of the system. However, by keeping the billiard area, the 



III. LEVEL STATISTICS 




(3) 




(4) 
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energy interval and the height and width of the teeth fixed, we find also for A 3 a smooth and systematic behavior, 
that depends on g. 

The behavior of A 3 (L) (see Figs. 2(a) and (c)) corresponds to the one of Iq. Those systems with small genus 
numbers g show high values of Iq and accordingly values of A 3 (L), which are close to L/15 (dotted line). Systems 
with high genus numbers g on the other hand, show Jo-values close to /Signer an( ^ accorc [ m g]y A3-curves, which are 
also closer to the Wigner- limit (solid line). Similar to the case of Seba billiards, they come very close to the Wigner 
distribution for small L. A 3 (L) turns out to be a more sensitive measure than Jo in the case of large (/-values. While 
the values of Jo for systems with genus numbers g > 20 lie already so close to the Wigner limit that a succession 
between them can hardly be recognized, we still observe clear differences between the individual curves of A 3 in the 
case of larger level distances. However, also the A 3 data indicate a systematic change from Poisson-like towards 
Wigner-like behavior with increasing g. 

In the following, we concentrate on I . We calculate Iq for several selected systems with different heights, widths 
and numbers of the teeth in the energy range of E S [0,3], which corresponds to roughly 40000 energy levels. The 
calculations were made in intervals of energy AE = 0.05, each of them containing about 600 — 800 levels. Here, 
we applied Neumann boundary conditions in all cases. The behavior of Iq under Dirichlet boundary conditions is 
qualitatively similar. 

First, we discuss the case of random values of the parameters <i, b X: by and h (Tab. II). For the first system R\ 
the values were uniformly distributed between Ad and 15d, for the second system i? 2 between Ad and lOd. The areas 
of the systems are again kept roughly constant. In Fig. 3 (a), Iq is plotted versus the energy E for the systems R\ 
(dashed lines) and R2 (solid lines). For most energy values, we find a roughly constant value of Iq that lies close to 
the Wigner- value. Only for very low energy values, E e [0,0.2] (which corresponds to roughly 2100 levels), Iq shows 
deviations from the high energy behaviour, towards higher values. This is due to the finite resolution of the boundary 
roughness in the limit of long wavelengths, where the systems look more regular. So, for a random structure of the 
boundary roughness, the high energy limit is reached quickly and no deviations of the Iq values from the Wigner limit 
are seen beyond the first 2000 states. 

In Fig. 3(b-c), I is plotted versus E for two types of systems with periodic boundary roughness, whose parameters 
are presented in Tab. III. In the first group (Fig. 3(b)), all systems have the same widths a = b x = b y = Ad and 
numbers N x and N y of the "teeth" (and correspondingly the same genus number g = 87). The heights of the teeth 
are different, h — Ad for system B\ (solid lines), 8d for system B 2 (dotted lines) and 16c? for system £> 3 (dashed lines). 
We first see that the narrow teeth account for a large low-energy regime with Poisson-like behaviour. Additionaly, we 
observe sharp peaks of the Iq values at several energies, which become more pronounced with increasing h and are 
most probably due to the periodic structures of the geometry. We will discuss these peaks in the following sections. 
The second group (cf. Fig. 3(c)) consists of billiards with larger widths a — b x — by — 8. This group of billiards shows 
less and only small peaks in the Jo-values and the high energy regime is reached quite quickly. 

Therefore we found two types of systems, where the high energy limit is reached quickly: geometries with broad 
teeth and geometries with a random distribution of the teeth widths. In systems with very narrow and regular teeth, 
on the other hand, energy windows appear, where the behavior of the system deviates significantly from the described 
behavior in the high-energy regime. In these intervals, the values of Iq are considerably larger than the expected 
high-energy values. These energy windows are most interesting and we look at them in more detail now. 

There are three different effects, which may lead to higher values of Iq . First, at smaller energies (larger wavelengths) , 
the teeth could not be sufficiently resolved and the states see only the rectangular main body of the system. This is 
the case for small energies, when half a wavelength is larger than the width a and the eigenfunctions are not small 
enough to penetrate the teeth. Second, also at higher energy values, the periodic structure of the teeth could allow 
for very regular functions, even if the boundary roughness is fully resolved. Also in this case, we expect a distribution 
close to Poisson. Third, also localized states can be a reason for the spectrum to behave Poisson-like. Therefore, we 
expect that Iq is closely related to special system properties and should be reflected in the shape of the eigenfunctions. 
In order to understand this, we now investigate the eigenfunctions in the different energy windows. 

IV. EIGENFUNCTIONS 

Some typical eigenfunctions of the system B3 under Neumann boundary conditions are presented on Fig. 4. They 
are taken from the different energy regimes (labeled from I to IV), as indicated in Fig. 4(f). Due to technical 
restrictions of the Lanczos algorithm, we could not calculate eigenfunctions of arbitrary high energies in large system 
sizes. The reason is that the density of states increases with the energy and the eigenfunctions become too close to 
each other and thus could not be seperated. At lower energies, we find the following characteristic shapes: 
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The eigenfunctions vf^ 7 ") and i]/( 7b ) of Fig. 4(a) and (b) are taken from the energy regime I, where Iq has the 
peak value of la = 0.942, very close to the Poisson value. The functions look different. While the function \f( 7a ) on 
Fig. 4(a) looks very regular and extended, function vf^ 76 ) on Fig. 4(b) is a rather localized (and regular) function, 
where non-zero amplitudes exist basically close to the boundary, i.e. inside the boundary teeth. Figure 4(c) (regime II 
on the Fig. 4(f)) represents rather the case of a chaotic function, i.e. the amplitude looks very random. Accordingly, 
we find Iq — 0.688 in this energy regime, closer to the Wigncr value. The function vI/( 777 ) on Fig. 4(d) (regime III on 
the Fig. 4(f)) looks again regular and Io = 0.893 is again close to the Poisson value. The difference of this function to 
the function \I>( 7a ) is the smaller wavelength in regime III. Here, the boundary roughness is to some extent resolved 
and the reason for the regular wavefunction lies in the periodicity of the teeth. These regular wavefunctions at higher 
energies should disappear in systems with random boundary roughness. Figure 4(e) (regime IV on the Fig. 4(f)) 
represents an interesting intermediate case of a function that is chaotic as well as localized. Iq in this case is 0.667, 
rather close to the Wigner value. The function \I>( 7V ') looks random in the inner rectangular part of the billiard, but 
its amplitudes in this part are very small. The largest amplitudes of the function are localized on the border of the 
teeth. So this case is in some sense intermediate between a localized and a chaotic function. This seems to be a weak 
localization mechanism, where the wave is reflected at the boundary roughness and interferes constructively inside 
the teeth. 

Accordingly, we found four characteristic types of eigenfunctions, regular ones, localized ones, chaotic ones and 
intermediate ones between chaotic and localized. All eigenfunctions seem to correspond to energy windows, which 
can be characterized by their corresponding Jo-values. We now introduce several measures for the eigenfunctions that 
allow us to distinguish between those cases. 

(i) As a first characteristics for the eigenfunction analysis we used the distribution of the amplitudes P(ip). By 
semiclassical arguments it was conjectured that for classically chaotic systems most eigenfunctions are a random 
superposition of plane waves, which leads to an amplitude distribution that is a Gaussian function [26-28], 

^) = ie^ (5) 



where ip — \^A^/ is normalized according to J \ip(x, y)\ 2 dxdy = A with the area A of the billiard. This normalization 
allows us to compare eigenfunctions of systems with different sizes. The amplitude distribution of eigenfunctions of a 
rectangular billiard, on the other hand, was shown to be [12]: 

P W = {^V)^(™)' °<M^ (6) 

where K(k) is the complete elliptic integral of the first kind. P(i/j) was already studied for the 7r/3-rhombus billiard 
with 5 = 2 [10] and for systems with a pointlike scatterer [12]. 

In Fig. 5 we show the amplitude distributions of our eigenfunctions from Fig. 4 by the open circles. (The filled 
circles will be explained in the next section.) The limiting cases of the Gaussian distribution for random functions 
and of Eq. (6) for regular functions are indicated by a dotted and a solid line, respectively. For the regular looking 
functions *< 7a ) and *< 777 ) (cf. Figs. 5(a) and (d)) we find a very good agreement with Eq. (6). The amplitude 
distribution on Fig. 5(b) for function \I/( 7b ) on the other hand consists of one large peak at ip s=s 0. This function 
is localized and only in a very small region of the billiard the amplitude is large. The function \f( 77 ) on Fig. 5(c) 
represents a function with a Gaussian distribution of the amplitude (cf. Eq. (5)). This is in line with our estimation 
that the function looks chaotic. The last case of the eigenfunction \|/( 7y ) that looks intermediate shows a distribution, 
which lies between the curves of the Gaussian and the localized function (cf. Fig. 5(e)). 

(ii) A second quantity to characterize the eigenfunctions is the localization volume (participation ratio) [29], 



r ( a ) = V 

A A J dxdy 



v ioc a — . ,,4 . . ' y > 



where is normalized according to J \^(x, y)\ 2 dxdy = 1 and V a — ^ / | SP^"^ | 4 dxdy^j . For some specific exampl 
of ^{x, y), we find 



Vioc = < 



1 for constant functions 

4/9 for regular sine or cosine-functions , . 

1/3 for Gaussian functions ^ ' 



k l/A for (^-functions. 



5 



For the functions of Figs. 4(a-e), we find V£ c a) w 0.429, V^> w 0.025, 3 w 0.332, w 0.441 and 

V£ V) w 0.060. V^ o) as well as v£* JI) is very close to the value of 4/9 of a regular cosine-function, whereas V loc 
is very small and corresponds to a localized state. The localization volume is very close to 1/3 for Gaussian 

functions, which confirms that the function is chaotic. For the last function, Vj Q ^ is again small, which means that 
also this function is rather localized, even if its V\ oc is slightly larger than In all cases, the values of V^J match 

very well to the amplitude distributions P(ip) of Fig. 5. 

(iii) As a third measure for the eigenfunctions, we investigate the energy surface. In order to do so, we expand their 
amplitudes C^m = (^ a; |$ n ,m) m the basis n,m of a rectangular billiard, which has the same linear extensions x,y 
as our rough billiard. Here, a enumerates the eigenstates of the rough billiard, while 3> n m are the eigenfunctions of 
the rectangular system. A similar analysis for chaotic billiards has been performed in Refs. [17,30]. 

The amplitudes |Cn™m| of our eigenfunctions are shown in Fig. 6. Except for the localized state ty (Ib \ where the 
functions & n .m(x,y) do not form a good basis, the values of |Cn™m| appear as peaks that are situated very close to 
the line of constant energy, E a <~ (n 2 + m 2 ). This means that those eigenstates $„ jrn (x, y) of the rectangular system, 
which have the energy E n ^ m close to E a , interfere and form the eigenstates of the rough system. However, the 
number of the participating states is very different. The regular states are represented in n, m-space by one large 
peak (one coefficient C^j has the absolut value close to 1), while the contribution of the others is vanishing. This is 

the case for the functions <]/( /a ) and \I»( 777 ), whose n, m-space can be seen in Figs. 6(a) and (c), respectively. \Cn,m\ 
of the chaotic state <]/( 77 ) is shown in Fig. 6(b). 

Here, we find many peaks of roughly equal heights along the energy surface, which shows that this state is a 
superposition of many different states m (x, y) that contribute roughly by equal weight and is therefore spread 
over the whole energy surface in n, m-space. A more complicated intermediate situation can be found in Fig. 6(d) 
(corresponding to the amplitude of ^ 7y )). Here, we find a dense distribution of peaks on the energy surface (as for 
the chaotic state), but also several secondary broader peaks that are not lying close to the curve of E a ~ (n 2 + m 2 ) 
and seem to be due to the high amplitudes localized in the "teeth", which cannot be described in the basis of 
the rectangular billiard. This state is obviously intermediate between localized and chaotic. Accordingly, also the 
n, m-space matches very well to the preceeding measures of the eigenfunctions and to the values of Iq. 



V. DISTRIBUTION OF EIGENFUNCTIONS 



We have seen in the last section that for five selected functions, all measures of the eigenfunctions match very 
well with the /o-results found from level statistics. In intervals, where the level statistics is closer to Poisson-likc 
behavior (I close to 1), we found eigenfunctions, which are either regular or localized. In intervals where the level 
statistics is close to Wigncr-like behavior (Iq w 0.637) on the other hand, the eigenfunctions seem to be chaotic with 
a Gaussian distribution of the amplitudes and spread over the whole energy surface in n, m-space. Now, we look at 
the distribution of the eigenfunctions inside a given energy interval. For this purpose, we calculated more than 100 
eigenstates for each of the four energy regimes of Fig. 4(f) (341 eigenfunctions in regime I, 131 states in regime II, 
139 states in regime III and 144 in regime IV). In each case, we calculated V\ oc and the amplitude distribution P(tp) 
over all eigenfunctions. 

In Fig. 7, we show the normalized histograms of the Vj c-values for the four energy regimes. In the first regime (see 
Fig. 7 (a)), we find two peaks: one peak is close to V\ oc = 0.44, which corresponds to eigenstates of a rectangular 
billiard, and one peak at V\ oc < 0.1, indicating localized states. In the second energy regime (see Fig. 7 (b)), we 
find a narrow distribution of V\ oc around the value of 0.33 of chaotic functions, in good agreement with the relatively 
low value of Iq = 0.688 in this regime. In the third regime (see Fig. 7 (c)), the values of V\ oc are distributed in the 
interval [0.35 — 0.46] with two peaks around 0.37 and 0.44. That indicates basically regular states with a slight end 
towards random behavior, which is in line with the value of Iq — 0.893 in this regime. In the fourth energy regime 
(see Fig. 7 (d)), most of the Vi oc -values are roughly distributed in the interval [0.1,0.2], which means again that the 
states are quite localized. Nevertheless the value of Iq = 0.667 in this regime is rather close to /Signer ^ ^^fa 
functions. This contradiction accounts for intermediate states between localized and chaotic. 

Next, we look at the average amplitude distributions in our four regimes, which are shown in Fig. 5 (a-e) by the 
full circles. Since in the regime I two types of eigenfunctions exist, we split them into two groups according to their 
Vioc-values. The average amplitude distribution for 263 states with Vi oc > 0.2 is presented in Fig. 5(a) and the one for 
78 functions with V\ oc < 0.2 on Fig. 5(b). The distributions in the regimes II — IV are shown in Fig. 5 (c-e) by full 
circles. The distributions are very close to those for the individual functions (open circles), from the preceeding section, 
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showing that the features, which we found for the single functions are characteristic for a whole energy regime. The 
eigenfunctions of regime I are either localized or regular, the eigcnfunctions of regime II chaotic, the eigenfunctions 
of regime 77/ again close to regular functions and in regime IV, there are a lot of intermediate functions between 
localized and chaotic. It is remarkable that even the intermediate functions cover a whole energy window and explains 
the discrepancy between the rather small values of V\ oc and the value of lo, which is close to the Wigner- limit. 

With these measures in mind, we can finally look at the systems with a random distribution of teeth widths, 
where lo approached the Wigner limit very fast. In Fig. 8, we show (a) the histogram of the Vi oc -values and in (b) 
the amplitude distribution of 152 eigenfunctions of system Ri in the energy interval [0.225 — 0.240]. The values of 
the localization volumes show a narrow distribution around the value of Vq oc = 0.33 of Gaussian functions and the 
amplitude distribution P(ip) coincides with a Gaussian function. So, these functions are quite obviously a random 
superposition of plane waves. 

VI. CONCLUSIONS 

In summary, we investigated the level statistics and the eigenfunctions of pseudointcgrable rough billiards with 
high genus numbers g. They had a rectangular body and a rough boundary with small "teeth" of different widths, 
distances and heights. The level statistics was found to be intermediate between Poisson- and Wigner-like behavior 
and approached the latter with increasing g. Especially for small level distances, the A3-curves can come very close 
to the Wigner statistics. This behavior was similar to the one of Seba billiards with an increasing number of pointlikc 
scatterers [12,14]. 

Additionally, we found different energy intervals with characteristic types of eigenfunctions, which correspond to 
a special behavior of the level statistics. In order to classify the eigenfunctions we employed several measures: the 
localization volume V\ oc , the amplitude distribution P{ip) and the behavior of the eigenfunctions in n,m-space. 

We found that all systems have a low-energy regime, where the wavelengths of the eigenfunctions are too large to 
resolve the boundary teeth and thus only see the main rectangular body of the billiard. The eigenfunctions in the low- 
energy regime have the characteristics of regular functions with the level statistics close to the Poisson distribution. 
However, for systems with a random distribution of the boundary teeth or for very broad teeth, this regime stays very 
small and the high-energy regime is quickly reached. 

For higher energies, the eigenfunctions of systems with cither random boundary roughness or with not too small 
teeth, are characterized by a close-to- Gaussian amplitude distribution P(ip) with a localization volume of V\ oc ss 0.33. 
They can be constructed by a random superposition of many regular functions. The second half moments Iq of the 
energy spacing distribution p(s) at higher energies are - apart from fluctuations - energy-independent and close to a 
Wigner-like behavior for large g > 20. Also the spectral rigidity A3 changes from Poisson- like towards Wigner-like 
behavior when g increases. However, A3 changes slower and still shows deviations from the Wigner behavior even for 
large g. 

For the case of periodic and narrow boundary teeth, we find many energy regimes, where the level statistics is 
close to Poisson behavior. We find that those functions are either localized with very small values of V\ oc or regular 
cosine functions with V\ oc close to 0.44. Also an interesting regime occurs, where the eigenfunctions show intermediate 
behavior with chaotic traces in the main body of the system but very high amplitudes inside the teeth, indicating 
weak localization. 

Accordingly, deviations from a Wigner-like behavior of the level statistics are due to three different effects, (i) a 
poor resolution of the boundary roughness at very small energies, (ii) regular wavefunctions due to periodic boundary 
roughness and (iii) localization. However, for random boundary roughness and not too small energies, those effects 
are not pronounced and the eigenfunctions can mostly be described by a random superposition of plane waves, which 
is for pseudointcgrable systems remarkable. 

For further research it would be interesting to investigate the eigenfunctions in the high-energy regime, which might 
help to understand the A 3 data for large values of g. As A 3 measures the long-range correlations between eigenvalues 
it seems likely that even a very small amount of regular or localized eigenfunctions in an energy interval of otherwise 
chaotic functions can lead to these deviations. 
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Geometry 
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8 


8 


21 


18 


13 


12 


49 



TABLE I.: Table of the geometries used in Fig. 2. The parameters a, h, b x and b y refer to the ones given in Fig. 1. The values 
of the parameters are given in units of the lattice constant d. 



Geometry 


a, b x , by, h 


Nx 


Ny 


9 


Ri 


[4, 15] 


22 


17 


77 


R2 


[4, 10] 


29 


22 


101 



TABLE II.: Table of the series of the geometries used in Fig. 3. The parameters a, b x , b y , h have random values from the 
indicated intervals. The values of the parameters are given in units of the lattice constant d. 



Geometry 


a = b x = by 


h 


N x 


Ny 


9 


Si 


4 


4 


48 


40 


87 


B 2 


4 


8 


48 


40 


87 


B 3 


4 


16 


48 


40 


87 


Bi 


8 


4 


24 


20 


43 


B 5 


8 


8 


24 


20 


43 


B 6 


8 


16 


24 


20 


43 



TABLE III.: Table of the series of the geometries used in Fig. 3. The parameters a, h, b x and b y refer to the ones given in 
Fig. 1. In both series of the geometries the width and distance between the teeth remain the same while the height is changing. 
The values of the parameters are given in units of the lattice constant d. 
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FIG. 1.: Sketch of the considered pseudointegrable geometry. The parameter a is a width of the "teeth", h is their height, b x 



and by are the distances between them in x- 



and y- direct ions. N x and N y are the numbers of "teeth" in x- and {/-directions. 



The genus number g of this geometry is g = 1 + Gi, where Gi is the number of salient corners with angles of 3n/2. The beam 
splitting property of these corners is shown by the two arrows that indicate two different trajectories of classical particles. 
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FIG. 2.: (a) and (c): The second half moments 7o are plotted versus the dimensionless energy E for several systems, (b) and 
(d): Aa(L) from the fixed energy interval E £ [1.5, 2.0] is plotted versus L for the same systems as in (a) and (c). Figs, (a) and 
(b) refer to Dirichlet, (c) and (d) to Neumann boundary conditions. The symbols refer to the systems of Tab. I with increasing 
genus number g = 3 (full cicles), g — 5 (open circles), g = 9 (full squares), g — 19 (a;), g = 29 (open squares), g = 39 (+) 
and g = 49 (full diamonds). In (a) and (c) the systems of g = 19 and 39 are omitted for a better visibility. Iq, A3 and L are 
dimensionless. In all figures, the Poisson-limit is indicated by a dotted and the Wigner-limit by a solid line. 
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FIG. 3.: The second half moment Iq is plotted versus the dimensionless energy E for Neumann boundary conditions and for 
different geometries (cf. Fig. 1). The systems of (a) have random parameters a, b x , b y and h with different widths of the 
distributions (cf. Tab. II). The solid line corresponds to the geometry R\, the dotted line for geometry R?. Figure (b) shows 
the systems Bi, B2 and B3 (cf. Tab. Ill), which have constant values of a — b x — b y — 4 and an increasing height: h = 4 (solid 
line), h — 8 (dotted line), h — 16 (dashed line). Figure (c) shows the systems B 4 , B 5 and B 6 (cf. Tab. Ill) with parameters 
a — b x — by = 8 and the height h — 4 (solid line), h — 8 (dotted line), h — 16 (dashed line). 
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FIG. 4.: (a-e) Typical eigenfunctions from different energy windows for the geometry B3. The amplitudes are indicated by 
different gray levels. The white regions stand for positive amplitudes, the black ones for negative amplitudes. The neutral gray 
tone stands for nearly zero amplitude. The black contour line shows the border and does not correspond to any amplitudes. In 
(e), the largest amplitudes lie at the border and are hidden by the contour line. In (f), Io is plotted versus the dimensionless 
energy E and the regions, from where the functions (a)-(e) are taken, are indicated by circles. The functions (a) and (b) (vt' /a ) 
and ^>' /i '') are both taken from the region I. 
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I I , I I I , I I I I I , I , I , I , I , l_ 

-2 2 -2 2 



FIG. 5.: The dimensionless amplitude distribution of (i) the single eigenfunctions presented on Fig. 4 (open circles) and (ii) the 
averaged amplitude distribution over many eigenfunctions in the corresponding energy intervals as explained in section V (full 
circles). For single eigenfunctions the letters (a)-(e) correspond to *< /a ), \&( Ib \ <I»( /7 ), ^i III ) > <E(( /V ) respectively. The dotted 
line indicates the Gaussian distribution and the solid one the distribution of regular sin or cos- functions (Eq. 6) . Functions (a) 
and (d) are regular; (b) - localized; (c) - Gaussian (chaotic) and (e) - intermediate between chaotic and localized. 
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FIG. 6.: Structure of the energy surface of the eigenfunctions (a) * (/a) , (b) * (//) , (c) # (///) and (d) from Fig. 4. We 

show the absolute values of the amplitudes |Cn"m| on the n, m-lattice, where a stands for the individual functions labeled la, 
II, III and IV. 
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FIG. 7.: The normalized histograms of the VI oc -values are shown for energy regimes (a) /, (b) II, (c) III and (d) IV. For the 
calculations were used 341 states in (a), 131 states in (b), 139 states in (c) and 144 states in (d). 
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FIG. 8.: (a) The normalized histogram of the Vioc-values is shown for the 152 eigenstates of the system _Ri from the energy 
interval [0.225 — 0.240]. (b) The averaged amplitude distribution of the same eigenfunctions as in (a) is indicated by full circles. 
The distributions for a regular sin or cos-function as well as for the Gaussian distribution are indicated by a solid and a dotted 
line, respecitvely. 
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